The American Naturalist

Male and female animals often differ remarkably in their reproductive behaviours, typified by competing males and caring females. One influential explanation traces sex-role divergence to anisogamy: sex differences in gamete size. This explanation remains controversial, because existing models have not fully resolved how anisogamy leads to secondary sex-role divergence, and exceptions to typical patterns are empirically well-documented. Here, we develop mathematical models linking anisogamy to secondary sex-role evolution, showing that anisogamy drives sex-role divergence along contrasting evolutionary trajectories depending on ecological context. When the sexes differ only in gamete size, male-biased mating competition and offspring care evolve under scramble competition, whereas uncertain paternity shifts the population towards typical patterns. By contrast, anisogamy generally drives typical sex roles under contest competition, although sex-role reversal can emerge as an alternative stable state when care must be apportioned among offspring. Anisogamy can thus explain the striking diversity of sex differences across taxa.

Demographic stochasticity is generally expected to contribute little to population variability when census abundance is large, because random demographic events are increasingly averaged across many individuals. This expectation, however, depends on how evenly individuals contribute to reproduction. In many natural populations, reproductive success is highly skewed, so that a large census population can coexist with a much smaller effective number of breeders (Nb). Although Nb is commonly used as a genetic measure of contemporary effective size, it also summarizes variation in reproductive contribution and may therefore contain information about demographic variability. Here, we develop a female-based age-structured theoretical framework linking Nb to temporal variability in mature population abundance, with N denoting the census number of mature females and Nb the effective number of female breeders. Recruitment variability arises from non-Poisson reproductive variation and environmental stochasticity, and its transmission to mature abundance is described by a life-history coefficient, {gamma}, determined by survival, maturation, and the contribution of overlapping cohorts. When the Poisson component of recruitment variance is negligible and Nb/N is sufficiently small, the coefficient of variation in mature abundance, CV[N], can be approximated from Nb, environmental variability, and {gamma}, with little additional information from census abundance. This approximation does not require specification of the functional form of the relationship between adult abundance and recruitment. Stochastic age-structured simulations closely reproduced the analytical prediction, even when additional demographic stochasticity in survival was included. Substantial fluctuations in mature abundance occurred at census sizes on the order of 105 when Nb was small, whereas overlapping cohorts buffered recruitment fluctuations. Because Nb can be estimated genetically from samples of a single cohort, the framework provides a direct link between genetic monitoring and population dynamics. More broadly, these results show that large census abundance alone does not guarantee weak demographic variability and highlight the potential of genetic data to inform expected population fluctuations when detailed demographic information is limited.

Understanding how species interactions, ecological drift, and immigration jointly shape community composition is a central goal of community ecology. These processes determine not only average community composition, but also how often different community states occur. Here, we ask how competition reshapes the neutral distribution of community composition generated by drift and immigration. We analyze a stochastic two-species competition model with immigration and derive the stationary distribution of relative species composition near neutrality. The stationary distribution shows when communities are usually balanced and when they are instead dominated by one species. Species dominance can become common even when the corresponding deterministic model predicts stable coexistence, particularly in smaller communities. Balanced and dominance states can also occur together, producing a trimodal distribution. Competitive asymmetry further biases the distribution toward the competitively favored species. Thus, the interplay between competition, ecological drift, and immigration determines not only how much community composition varies, but also which community states are most likely to be observed. Our results provide an analytical bridge between neutral and competition theory and a framework for interpreting variation in community composition.

Haldane's principle relates equilibrium mutation load to deleterious mutational input rather than to selection strength. We study this balance in a homeostatically regulated, damage-structured population where mutation occurs during reproduction, repair restores individuals toward lower-damage states, and nonlinear demographic feedback maintains constant population size. An exact demographic identity shows that load at every time equals excess reproductive effort weighted by mean relative fecundity. Without repair, the equilibrium branch with positive pristine-class frequency is exactly solvable, and along it the load is exactly the mutation probability times the realized per-capita reproductive turnover rate; when the terminal cap is removed, this is precisely the realized deleterious input, for every mutation probability below one. For rare mutation it reduces to the mutation probability times the pristine turnover rate, with the local fitness contrast canceling at leading order. With a finite terminal class the branch instead reaches a closure-dependent boundary collision. Repair modifies the rare-mutation balance by partitioning first-order removal from the first damaged class between selection and restoration, multiplying the leading load by the selective share of that removal. An exact binary solution extends these results to finite mutation pressure: the unrepaired interior equilibrium collides nonhyperbolically with the fully damaged boundary and converges algebraically at the threshold, whereas any positive repair rate restores a unique interior equilibrium for positive mutation; at the no-repair threshold, its displacement from the boundary scales as the square root of the repair rate as repair vanishes. These results separate the familiar competing-transition structure of mutation--selection balance from the specifically homeostatic effects of turnover-dependent mutation input and nonlinear demographic feedback.

Assortative mating for a trait under selection can increase the rate of adaptive evolution and limit introgression of detrimental alleles. We therefore hypothesized that assortative mating facilitated the rapid evolution of flowering phenology observed in the self-incompatible perennial Lythrum salicaria (purple loosestrife) as it invaded eastern North America. We characterized flowering schedules of 600 individuals representing 20 populations sampled along a latitudinal gradient of 10 degrees and grown in a common garden at mid-latitude. Consistent with our hypothesis, we found that flowering schedules were less synchronous in the mixed population of the common garden (s = 0.50) relative to the 20 locally adapted populations comprising it (s = 0.56 to 0.86), and flowering synchrony decreased towards higher latitudes of origin. Despite high synchrony, assortment potential ({rho}) by growth and phenology varied by population ({rho} = 0.03 to 0.61), but the predicted negative correlation between s and {rho} was masked by among-individual fitness variation in the common garden. We show how fitness variation and asynchrony associated with locally adapted phenology can reduce the fitness impacts of maladaptive gene flow and accelerate the rate of adaptive evolution during invasion.

Ecological stability and function emerge from the many and diverse types of interactions among species. Ecological networks represent the structure and strength of species interactions, but their study is often restricted to single interaction-types, likely due to an absence of models for plausibly (i.e., nonrandomly) structured 'multiplex networks' of multiple interaction types. Here, we present a general but mechanistic model that simulates multiplex networks using the different, but correlated, traits that determine different species interactions. We show that ecological stability varies systematically with these trait correlations by controlling the emergent complexity and positioning of strong interactions in the community. Our model provides a flexible but power tool to study multiplex networks under alternative trait-based hypotheses, offering new insights into the structure and stability of complex natural systems.

Shifts in predation and human harvest often coincide with changes in life history trait distributions in prey populations. These changes may arise via evolution or plasticity, and extrinsic environmental effects may affect such traits. This complexity prompts questions about how trait distributions are shaped by multiple factors' interactions. We hypothesized that food availability and climate affect both the magnitude, and mechanism (evolution vs. plasticity), by which traits change in response to harvest's intensity and size selectivity. We tested these hypotheses using individual-based bioenergetic population models allowing both plasticity and evolution. We exposed genetically diverse fish populations to gradients of productivity and stochasticity in food availability, crossed with gradients of selective or nonselective fishing mortality. Food availability shaped traits both plastically and evolutionarily, and modulated how harvesting influenced those responses. Food effects can suppress harvest impacts altogether, or obscure evolution's effects on phenotypes. However, when climatic factors shape food-driven changes by affecting both fish and their prey, either synergistic or antagonistic interactions can emerge between food variation and harvest in determining phenotype distributions. Our research suggests that understanding consumer-resource dynamics, and climatic processes shaping them, is essential to interpreting how harvest can shape traits relevant to population viability and fisheries sustainability.

Defended prey are often thought to rely on conspicuous warning signals to distinguish themselves from undefended prey, facilitating predator learning and avoidance. Yet many chemically defended organisms are inconspicuous and difficult to distinguish from palatable prey, raising a fundamental question: how do predators avoid defended prey when prey themselves provide little or no reliable signals of defense? We tested whether host plants can function as indirect indicators of prey defense using Carolina wrens (Thryothorus ludovicianus), chemically defended Eumaeus atala larvae, undefended Galleria mellonella larvae, and two distinct plant species. Defended and undefended prey were experimentally made to appear visually identical, forcing predators to rely on host-plant associations as the only available source of information. Wrens rapidly altered both search behavior and attack decisions according to plant identity. Mechanistic attack models incorporating predator discrimination between defended and undefended prey were strongly favored when plants were present, with defended prey experiencing substantially lower attack risk than undefended prey. When plants were removed, support for discrimination weakened and attacks on both prey types increased. Spatial autocorrelation of host plants further shaped predator behavior, with strongly clustered plants hosting defended prey receiving disproportionately fewer visits. Together, these results suggest host plants can function as ecologically meaningful cues of prey defense, influencing both predator search patterns and attack decisions. More broadly, our findings suggest warning systems may extend beyond prey phenotype alone, and that ecological context (including host-plant identity and spatial structure) may contribute to the evolution and maintenance of aposematism.

Environmental variability influences the evolution of phenotypic plasticity and the ecological interactions among species in communities, but it is unclear how these processes interact to affect ecological systems across scales. We investigate models describing how plasticity evolves in competitive communities exposed to seasonal variation in abiotic conditions. We conducted eco-evolutionary analyses to identify communities of coexisting species with evolutionarily stable plasticity traits. These analyses revealed that the presence of interspecific competitors constrained the evolution of plasticity in some species, which instead evolved to be specialists. Plasticity often impaired species coexistence by allowing individual species to control more niche space. Communities containing plastic species were usually less species-rich but more productive and less variable. Similar results were obtained when plasticity rate was evolutionarily labile, but plasticity rate altered the patterning of temporal niche partitioning. This work shows how plasticity and competition combine to shape ecological and trait responses to environmental variability.

Functional traits are usually assigned to individuals, yet some properties exist only as differences among population members. Whether such higher-order properties predict collective performance remains unclear. Here we used the aquatic plant Spirodela polyrhiza to test whether phenotypic and genetic differences among genotypes predict non-additive population growth. We grew twelve Japanese strains alone and in all pairwise mixtures and quantified their growth dynamics. Mixtures sometimes outperformed the mean of their component monocultures, and some exceeded the better-performing monoculture, although diversity effects were negative on average. Positive effects were most strongly associated with between-strain differences in frond colour and morphology, identifying specific dimensions of phenotypic diversity as population-level predictors of synergistic growth. Genome-wide higher-level association mapping identified loci where between-strain diversity was associated with growth effects, implicating stress responses, secondary metabolism and frond morphology. Loci associated with reduced performance implicated chemical perception and self/non-self recognition, whereas those associated with positive effects reflected functional differentiation in stress responses and physiological niches. Our findings suggest that some traits generate new functions at the population level through variation among individuals, making diversity itself a higher-order functional property.

Evolutionary rescue occurs when genetic adaptation prevents extinction of a declining population. Mutation and recombination generate adaptive variation, but when their rates at adaptive loci are too low, rescue may require evolution toward higher rates of genetic change. Here, we derive and compare analytical approximations for the probability of rescue through two alternative routes: mutation-rate and recombination-rate evolution. Our model comprises two resident genotypes from which the same adaptive haplotype can arise along either route, with rescue requiring the appearance of an unlinked modifier allele that increases mutation at, or recombination between, two fitness-determining loci. We show that resident genotype composition determines the relative success of the two routes. Mutation-rate evolution yields a higher rescue probability when one resident genotype is sufficiently more common because mutation can act on either genotype whereas recombination requires both. Otherwise, the two routes yield similar rescue probabilities. Even when either route makes rescue appreciably likely, however, the modifier alleles facilitating rescue are often lost before rescue is complete. This loss is driven by multiple factors: segregation can separate the modifier from the adaptive haplotype; the fitness cost of modifier expression raises extinction risk for modifier-bearing lineages; and continued modifier-induced mutation disrupts the adaptive haplotype after it forms. The latter effect can make recombination-rate modifiers more likely to persist even when mutation-rate evolution yields a higher rescue probability. Overall, these results show that both mutation- and recombination-rate evolution can induce rescue, but successful rescue need not entail persistence of the modifier allele that facilitated it.

Quantifying how individuals adjust the seasonal timings of life-history events (i.e. their phenology) in response to seasonal environmental conditions is fundamental to predicting impacts of changing seasonality. However, the current standard approach of inferring phenological responses by directly analysing observed event dates relies on foundational assumptions that could be invalidated by changing seasonality, and excludes individuals that do not express focal events, potentially generating misleading inferences. To overcome these limitations, we devised a Phenological Threshold Model (PTM) that conceptualizes phenological events as switches between discrete life-history states and primarily estimates variation in underlying latent dynamic responses to seasonal progression, from which event dates secondarily arise. We demonstrate that the PTM can accurately retrieve simulated parameter values, including given incomplete seasonal field observations. Application to multi-year data on annual breeding in black-legged kittiwakes (Rissa tridactyla) demonstrated substantial among-individual variation in latent responses to seasonally varying environmental drivers, and associated daily probabilities that individuals will initiate breeding, revealing how annual distributions of breeding dates and proportions of non-breeders jointly arise. We thereby formalise a new conceptual and statistical framework for quantifying the potential for, and constraints on, rapid plastic and micro-evolutionary changes in phenology in changing seasonal climatic and environmental regimes.

Offspring can inherit a physiological or molecular state from their parents, and the carry-over in phenotypes across generations that results is often called epigenetic or phenotypic "memory''. Most theory on the evolution of such memory asks how the rate of switching between pre-existing phenotypic variants evolves in a fluctuating environment. Here we ask instead how the variants themselves originate. We model an asexual population in a periodically alternating environment in which a genotype specifies two inherited states, the phenotypic value each produces and the probability with which each is retained from parent to offspring. Starting from an ancestor that produces a single phenotype, mutations creating small differences between inherited states have no first-order selective advantage. Their spread depends on second-order selection and requires strong selection against phenotype-environment mismatch. A maternal effect of the same small size, conditioned on the maternal environment, is under first-order selection in the same model. Once phenotypes have diverged, genotypes whose retention probabilities match the duration of the environmental states are strongly favored, and the evolved strategy persists under conditions in which it could not originate. Memory is a consequence of individual selection and its origin probably follows earlier selection for cue-dependent maternal effects.

Ecological uniqueness, commonly quantified by local contribution to beta diversity (LCBD), is important for understanding and conserving biodiversity. Its widely documented negative relationship with species richness has led to the view that species-poor communities may be disproportionately important for regional diversity. Here, we derive neutral expectations for richness-uniqueness relationships and show that negative slopes can arise from neutral assembly rather than necessarily reflecting an ecological signal of uniqueness. Neutral simulations validate the analytical expectations across abundance- and incidence-based LCBD formulations. Applying the framework to the Barro Colorado Island tropical forest and North American Breeding Bird Survey, we show that accounting for neutral sampling expectations substantially changes the ranking and spatial distribution of ecologically unique sites. By separating neutral sampling expectations from excess uniqueness, our framework provides a principled basis for interpreting ecological uniqueness and informing biodiversity conservation.

Extreme climatic events (ECEs) could potentially induce episodes of rapid adaptive evolution, altering the evolutionary trajectories and viability of wild populations experiencing increasingly stochastic environments. Adaptive potential stems substantially from additive genetic variance (V) in survival. However, the degree to which ECEs impact this variance, and the resulting evolutionary implications, have not been fully conceptualised or quantified in nature. These aims are highly challenging because underlying genetic and environmental variances have fundamentally non-linear effects on the dichotomous phenotypic outcome of survival versus death, while substantial data are required for estimation. We present a conceptual and analytical framework that formalises biological hypotheses for how V depends on environmental severity, how such dependences are reshaped by non-linear genotype-environment-phenotype relationships and demographic structure, and the implications for evolutionary responses. We then use 25 years of survival and pedigree data from European shags (Gulosus aristotelis) to decompose variance in adult survival across diverse environmental conditions including ECEs. Using a capture-recapture animal model, we demonstrate non-negligible V underlying adult survival across all years combined, implying notable adaptive potential. Yet, partitioning by environmental severity revealed substantial V in good years with high survival but negligible V in bad years with ECE-induced low survival, implying a genotype-by-environment interaction that constrains adaptive potential under severe conditions. However, scale transformations showed how the non-linear genotype-environment-phenotype relationship can distort or even reverse such patterns, by reducing the extent to which V is phenotypically expressed in good conditions with high mean survival and amplifying it in more severe conditions with lower mean survival. Evolutionary responses to severe conditions could consequently be non-negligible even given small underlying V. We thereby demonstrate how ECEs can reshape adaptive potential, highlighting the need to integrate impacts of environmental severity on genetic variation in fitness components into predictions of population resilience under climate change.

Beta diversity is a fundamental component of regional biodiversity that accumulates from compositional differences between local species assemblages. Although beta diversity is increasingly applied across spatial scales, its scale dependence remains incompletely synthesized, particularly concerning variation in spatial scale used to define alpha diversity (grain size). Using tree community data from Barro Colorado Island, we find that both additive and multiplicative beta diversity are consistently higher for unequal than equal grain divisions across different spatial extents and numbers of grains. These results can be explained by Jensen's inequality based on the concave-down relationship of alpha diversity with area that in turn increases beta diversity through the scaling relationship between alpha and gamma diversity. Combined with our additional analysis showing that this prediction is robust to more general models of species richness, including species--area and species--sampling-effort curves, these results reveal a previously unrecognized, predictable range of beta diversity variation.

Animals encounter a variety of food items and must determine their future value and preservability. Because animals in laboratory settings are typically provided with uniform food items, there is a dearth of knowledge on how animals assess natural complex food. Here, we dissect how wild agoutis-a critical seed disperser of the Neotropics-decide whether to consume or bury food in small caches, i.e. scatter-hoarding. We find that agoutis preferentially cache nuts over fresh fruits and, when foraging, undergo a sharp transition from eating to scatter-hoarding. Remarkably, agoutis scatter-hoarded empty nuts with intact shells but not filled nuts with broken shells. This decision was dependent on an oro-manual food assessment, an ecological information-seeking behavior. They used their lips and incisors - a somatosensory fovea - to navigate the surface of the nut optimally accumulating decision-critical sensory information. Together, these findings provide a mechanistic account of preservability assessment in a wild rodent, revealing how information-seeking behaviors support long-term value optimization in natural environments.

The Yampolsky-Stoltzfus (YS) model provided the first formal demonstration that a bias in the introduction of new mutations can shape the direction of adaptive evolution, even when mutation rates are small relative to selection coefficients. Despite its foundational role, the model has been characterized almost entirely by simulation and closed-form solutions for its central quantity, the probability PA that a mutationally favored but less-fit genotype fixes ahead of a rarer, fitter alternative, have been available only in the limiting regimes of vanishing mutation supply and infinite population size. Here we derive closed-form approximations for PA across the intermediate regime of finite population size and appreciable mutation supply that lies between these two limits. Under strong selection between the competing mutants, we combine an establishment-probability argument with clonal interference theory to obtain an explicit expression for PA that accounts for both the ongoing input of the mutationally favored mutation and the frequency-dependent establishment probability of its competitor. Under weak selection between competing mutants, we derive a complementary expression based on the fixation of a mixed population under Kimura's classical formula. Both expressions agree closely with direct simulation of the underlying Wright-Fisher process, recover the known limiting cases, and remain consistent with one another over a substantial intermediate range of parameters. These results are expected to provide a basis for assessing how the strength of mutation bias affects adaptation in empirical data.

In natural ecosystems, community dynamics are driven by diverse interactions operating through both ecological and reproductive processes. Although theory predicts that these processes jointly generate complex and contingent dynamics through the interplay of stabilizing (negative frequency dependence arising from niche differences) and destabilizing effects (positive frequency dependence caused by reproductive interference), empirical tests of this prediction remain scarce. Here, we combined experimental and theoretical approaches to investigate the dynamics of two congeneric annual plant species (Commelina) experiencing both ecological competition and reproductive interference. Our empirically parameterized model revealed a large niche difference and mutual reproductive interference. However, the reproductive interference was not strong enough to generate pronounced contingency of community dynamics on initial conditions. Instead, coexistence was the most likely outcome regardless of initial species composition. Parameter exploration indicated that the weak contingency arose because reproductive assurance strategies of Commelina species buffered the minority disadvantage caused by reproductive interference. These findings provide empirical evidence that reproductive interference does not necessarily destabilize coexistence and that it may be outweighed by ecological niche and fitness differences. Our study highlights the importance of integrating ecological and reproductive processes towards better prediction of community dynamics.

Flowering plants vary in sex allocation (resource division between male and female function) and in selfing rate. These traits should evolve together: selfing changes fitness returns from pollen and ovules, whereas sex allocation changes the pollen available for selfing and export. Yet most models allow only one trait to evolve, leaving their evolutionary feedback unresolved. Here we model the joint evolution of sex allocation, pollen export and maternal control of competition between self- and cross-pollen. We show that selfing and sex allocation strongly influence one another, as selection favours associations between sex allocation and mating-system traits. When maternal control is fixed to recover classical mass-action, evolution can lead to the stable coexistence of selfing hermaphrodites and males specialised in pollen export. This androdioecy persists because hermaphrodites cannot prevent fertilisation by males. When maternal control of pollen competition also evolves, we find that where a mixed-mating equilibrium exists, selection again favours the emergence of androdioecy. As hermaphrodites evolve to favour their own pollen, however, they progressively exclude pollen from males, driving them extinct. The population may then return to mixed mating and repeat the process, or evolve complete selfing. We discuss these results with reference to known variation in the mating and sexual systems of plants.